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THE ORIENTED SWAP PROCESS 

By Omer Angel/ Alexander Holroyd^ and Dan Romik 

University of Toronto, University of British Columbia and The Hebrew 

University 

Particles labelled 1, . . . ,n are initially arranged in increasing or- 
der. Subsequently, each pair of neighboring particles that is currently 
in increasing order swaps according to a Poisson process of rate 1. 
We analyze the asymptotic behavior of this process as n — > oo. We 
prove that the space-time trajectories of individual particles converge 
(when suitably scaled) to a certain family of random curves with two 
points of non-differentiability, and that the permutation matrix at a 
given time converges to a certain deterministic measure with abso- 
lutely continuous and singular parts. The absorbing state (where all 
particles are in decreasing order) is reached at time (2-f o(l))n. The 
finishing times of individual particles converge to deterministic lim- 
its, with fluctuations asymptotically governed by the Tracy- Widom 
distribution. 

1. Introduction. Let 5„ be the symmetric group of all permutations a = 
(o"(l), . . . , a{n)) on {1, . . . , n}, with composition denoted ((Tr)(z) := a{T{i)). 
For 1 < i < n — 1, denote the adjacent transposition or swap at location i by 
Ti := {i i + 1) = {1,2, ... ,i + l,i, ... ,n) G Sn- 

The oriented swap process is the 5„-valued continuous-time Markov pro- 
cess {r]t)t>o = {v?)t>o defined as follows. The initial state rjQ is the iden- 
tity permutation id := {1,2, ... ,n). From a state a, for each i satisfying 
a{i) < a{i + 1), the process jumps at rate 1 to the state crrj. Note that the 
reverse permutation rev := {n, ... ,2,1) is the unique absorbing state. Our 
focus is on the limiting behavior of the process (?/")t>o as n — > oo. We call 
the random permutation r]t the configuration at time t, and we call ri'j~^{k) 
the location of particle k at time t. We call the function t^rj'[^{k) the 
trajectory of particle k; see Figures 1-3. 
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Fig. 1. An illustration of an oriented swap process with n — 5. Trajectories are shown 
by lines. 

Our first result states that the trajectories converge to a certain family 
of random curves. The limiting curve for a particle at a given location is 
deterministic, once a random initial speed has been chosen, and is smooth, 
except at two points. Define the scaled trajectory TJ! = Tj^ : [0,oo) — > [0, 1] of 
particle k by 




500 1000 1500 2000 

(a) 




0.5 1 1.5 2 0.5 1 1.5 2 



(b) (c) 

Fig. 2. (a) Selected particle trajectories in a simulated oriented swap process with 
n = 1000; (b) selected possible limiting trajectories for particle [3n/10j; (c) selected lim- 
iting trajectories (see Theorem 1.1). 
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(b) 



Fig. 3. (a) The configuration (i.e., the support of the permutation matrix) at times 
0, n/3, . . . , 6n/3 for a simulated oriented swap process with n = 500; (b) an illustration 
of the limiting measures for these configurations as n —> oo (see Theorem 1.2). 



Theorem 1.1 (Trajectories). Let k = k{n) he a sequence satisfying k/n^ 
y € [0, 1] as n— >oo. Then the scaled trajectory of particle k satisfies 

Tfc =^ (j)y as n ^ oo. 

Here, denotes convergence in distribution with respect to the uniform 

topology on functions [0, oo) — > M, and (j)y is a random function given by 

L-{s)V {y + Us) AL+{s), s<-fy, 
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where U is uniformly distributed on [—1,1], and we have the deterministic 
functions 



L-{s) ■.= y + s-2^, L+{s):=y-s + 2Js{l-y) 



jy:=l + 2^y{l-y). 

Above, V and A denote maximum and minimum, respectively. Note that 
oV6Ac=(aV6)Ac = oV(&Ac) is well-defined when a<c and is the point 
in [a, c] closest to b. 

The configuration at a given time is a random permutation. The next 
result states that, after suitable scaling, its permutation matrix converges 
to a certain deterministic limit. Define the scaled configuration ^" = fig at 
(scaled) time s to be the following random measure on [0, 1]^: 



(1) 



1 



n 



k=l 



n n 



where 6{x, y) denotes the point measure at (x, y) E M^. The limiting measure 
has an absolutely continuous part, and a singular part which concentrates on 
a curve. For functions /i : M — > M, we denote by h{x)6y^gf^^^ dx the measure 
on which concentrates on the curve {{x,y):y = g{x)} and assigns measure 
h{x) dx to Ax R. 

Theorem 1.2 (Configurations). For any s > 0, the scaled configuration 
Us satisfies 

/i" =^ Ks as oo. 

Here, denotes convergence in distribution with respect to the weak 

topology for random Borel measures on R? , and Ks is the deterministic mea- 
sure on [0, 1]^ given by 

Ks{dx X dy) 



(2) 



+ 1 



L-{s)<y<L+{s] 
X — s<y<x + s 

s>l, 
x^{W-,W+) 



— dxdy 
2s ^ 



-'y=l—x 



dx 



s<l, 
X < s 



+ 1 



s > 1, 

W- <x< W: 



6 T-i \dx 



+ 1 



s < 1, 
X > 1 — s 



+ 1 



.S>1, 

Wr <x< W+ 



dx, 



where are as in Theorem 1.1, and for s > 1, we write Wf^ := i=*=v^ 
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The functions W and 7 from Theorems 1.1 and 1.2 are inverses; more 
precisely, we have 7^± = s for all s £ [1,2]. 

The limiting measure Kg in Theorem 1.2 is symmetric under the transfor- 
mations {x,y) I— > (1 — x, 1 — y) and {x,y) 1— > {y,x). The former corresponds 
to an obvious symmetry of the model, but the latter reflects the somewhat 
surprising fact that the permutation rjt is asymptotically symmetric in law 
under inversion. In fact, this symmetry holds exactly for each n and t, as 
the following result states. 

Theorem 1.3 (Symmetry). For any n and t, the permutation rj^ is 
equal in law to its inverse (ryf)"^. 

The inversion number of a permutation cr is defined as 



An alternative description of the process (r/j) is that the jump-rate from a 
to (TTj equals l[inv(crTj) > inv((T)]. From Theorem 1.2, we can deduce the 
following: 

Theorem 1.4 (Inversion number). For each s > 0, the scaled inversion 
number of the configuration satisfies the convergence in probability 



The limiting function on the right-hand side above is analytic, except at 
s = 1 and s = 2, where it is, respectively, three times and once continuously 
differentiable. 

Remarks on time-parameterizations. One may consider several natural 
variants of the process {r]t) in which time is parameterized differently. In the 
version introduced above, the total jump rate from a permutation a equals 
the size of its ascent set A = A{a) := {i:a{i) < a{i + 1)}; we refer to this 
as the variable- speed continuous-time process. In the fixed-speed continuous- 
time process, at rate 1, we choose an i uniformly from A{a) and jump to 
(7Tj. In the variable-speed discrete-time process, at each step, we choose a 
uniformly random i from {1, . . . ,n — 1} and jump to <TTj, provided i £ A; 
while in the fixed-speed discrete-time process, we choose i uniformly from 
A. Clearly, the sequence of distinct states visited has the same law for each of 
these four processes, and the ith state to be visited always has inversion num- 
ber i. Therefore, using Theorem 1.4, one may easily translate Theorems 1.1 
and 1.2 into analogous asymptotic results for the other three processes, with 



inv(cr) := #{{ij):i <j and a{i) > a{j)}. 
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the limiting objects being identical, except for a deterministic time change 
(for brevity, we omit the full statements). Surprisingly, however, the exact 
symmetry in Theorem 1.3 applies only to the variable-speed models. For ex- 
ample, for the fixed-speed discrete-time process, with n = 4 and at time step 
3, it is easy to check that the two mutually inverse permutations (2,4, 1,3) 
and (3,1,4,2) have respective probabilities 1/3 and 1/6. 

Lastly, define the finishing time /3"(/c) = (3{k) of particle k to be the 
(random) last time at which it moves, reaching its final position: 

(3{k) := sup{t > : 7]t{k) j^n + l-k}. 

Also, define (3^: = (3^ := maxk f3{k), that is, the hitting time of rev. The fol- 
lowing result is strongly suggested (although not directly implied) by The- 
orem 1.1: 



Theorem 1.5 (Finishing times). We have the convergence in probability 
max|/3"(A;)/n-7fc/„| ^0 



as n — > oo, where 7^ := 1 + 2y/y{l — y) is as in Theorem 1.1. In particular, 
we have 



/3:/nA2. 



Furthermore, we establish the following result on the fluctuations of the 
finishing times from their typical values. The Tracy-Widom distribution 
function is -Ftw(-2) := exp[— /^°°(x — z)u{x)'^ dx], where u{x) is the unique 
solution of the Painleve equation u" = 2u^ + xu having the Airy function 
asymptotics u{x) ~ 1/ {2y/TT)x^^/'^e~'^^^^^ as x ^ 00. The Tracy-Widom 
distribution originally arose [14] in random matrix theory, and has since 
been found to appear as a limiting law in several combinatorial models; see 
[3, 5, 8]. 

Theorem 1.6 (Finishing time fluctuations). Let k = k{n) be such that 
k/n^ y £ (0,1) as 00. We then have the convergence in distribution 

P''{k)-7k/nn 

— ^ -TTW as n—> 00. 



7,'/'(y(l-y))-V6nV3" 

Here, the requirement that y E (0, 1) is needed. For example, the finishing 
time /9"'(1) of the first particle has Gamma(n — 1,1) distribution (since the 
particle jumps only to the right, always at rate 1), which converges, after 
scaling, to a Gaussian limiting law. 



THE ORIENTED SWAP PROCESS 



7 




I 



Fig. 4. Selected trajectories and half-time configuration for a 2000-element uniform sort- 
ing network. 

Remarks on sorting networks. A sorting network is any sequence (ii, . . . , ZAf) 
of length := (2) = inv(rev) such that Tj^ • • • Tij^ = rev. Clearly, the sequence 
of swaps in the oriented swap process (r/f ) corresponds to a random element 
(with a certain non-uniform distribution) in the set of n-particle sorting 
networks. The uniform sorting network is instead chosen according to the 
uniform distribution on the same set. The present work was in part moti- 
vated by the striking results and conjectures on uniform sorting networks 
in [2]. The two processes behave quite differently, but share some features. 
For example, in the uniform sorting network, the particle trajectories con- 
jecturally converge to random elements in a one-parameter family of curves, 
while the configurations conjecturally converge to a family of determinis- 
tic measures (sine curves and projected sphere measures, respectively; see 
Figure 4). Both of these properties hold for the oriented swap process (The- 
orems 1.1 and 1.2 above), but with different limiting objects. 

Remarks on the proofs. Our analysis of the oriented swap process relies 
on a connection with the theory of the totally asymmetric simple exclusion 
process [TASEP). We show that the oriented swap process can be rep- 
resented in terms of a family of coupled TASEPs, by observing that the 
behavior of particles 1, . . . , fc (if we ignore their labels) is that of a TASEP 
on the finite interval [l,n] and then representing the TASEP on the finite 
interval in terms of the TASEP on Z using a combinatorial mapping. The 
main probabilistic results will follow using known limiting results for the 
TASEP. For Theorems 1.2 and 1.4, we use the classical hydrodynamic limit 
theorem of Rost [13] for the TASEP. Theorem 1.1 uses a result of Guiol and 
Mountford [11] on the trajectory of a second-class particle in the TASEP. 
Theorem 1.6 is a consequence of a theorem of Johansson [8] on the conver- 
gence of the scaled fluctuations of percolation times in oriented last-passage 
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Fig. 5. The permutahedron for n = 4. 

percolation in to the Tracy- Widom distribution. See also [1, 6] for results 
related to the joint distribution of trajectories. 

Further remarks. The permutahedron is a natural embedding in Eu- 
clidean space of the Cayley graph of Sn with the nearest-neighbor swaps 
Ti, . . . , Tn-i as generators. Specifically, the vertex corresponding to a permu- 
tation a is located at the point = {a~^{l), . . . , a~^{n)) G M"; it is easily 
verified that all vertices lie on an (n — 2)-sphere (see Figure 5). The oriented 
swap process may be regarded as a continuous-time simple random walk on 
the permutahedron, constrained to move only in directions which increase 
the scalar product with the vector from id~^ to rev~^. 

The oriented swap process is a special case in the family of Markov pro- 
cesses on Sn in which the jump rate from a to fXTj is a if cj(i) < a{i + 1), and 
a' otherwise, for fixed parameters a,a' . These processes were studied in [4], 
where it was proven that the mixing time for a 7^ a' is 0{n) [or O(n^) in the 
discrete-time version of the process of that paper]. Analogs of Theorems 1.2 
and 1.4 for these processes can be proven along the same lines as our proofs, 
by using results on the partially asymmetric exclusion process due to Kip- 
nis, 011a and Varadhan [9, 12] that generalize the results we used for the 
TASEP. Somewhat surprisingly, the symmetry condition in Theorem 1.3 
also holds for these models; see [1]. However, extending Theorem 1.1 to the 
partially asymmetric swapping process would require proving an analogue of 
the Ferrari-Kipnis theorem for partially asymmetric exclusion; see Section 
8. 

2. Symmetry. In this section, we prove Theorem 1.3. As noted, the result 
holds for both the discrete- and continuous-time variable-speed versions of 
the model. The continuous-time claim follows from the discrete-time state- 
ment, which, in turn, is an immediate consequence of the following lemma. 
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Define the operator Si, acting from the right on permutations, by 

^_Q_(crTi, if (T(i) < cr(i + 1), 
* ' \a, otherwise 

for fj G 5„; thus, Si sorts the particles in positions i and i + 1 into decreasing 
order. Note that we can write 

a • Si = max{cj, ctTj}, 

where the maximum of two permutations is the one with the greater inver- 
sion number. Note that inv(cjTj) — mv{a) = ±1, so the maximum above is 
always well-defined. 

Lemma 2.1. For any sequence we have 

(3) id- Si,- --Si^ = {id - Si^ - - - Si,y^ . 

Proof of Theorem 1.3. The discrete-time variable-speed process may 
be constructed in terms of an i.i.d. sequence of random variables {Ij)j>o cho- 
sen uniformly at random from {1, . . . ,n — 1}. The process is then given by 
Cfc = id • Sj, - - - . Since reversing Ii, . . . ,1^ does not change their joint dis- 
tribution, the theorem follows from Lemma 2.1. The continuous-time process 
may be defined as 77" = Cx(t)^ where X(-) is the counting function for a Pois- 
son process with rate n — 1 which is independent of the swap locations. The 
result follows. □ 

Proof of Lemma 2.1. We first prove the following claim. Let . . ,{„ 
be any integer sequence. If we let X = id • iSj^ • • • 5^,^ and Y = Ti - Si^^ - - - Si^ , 
then 

y-^ =max{X^^X^V^} = X-^ • Si. 

Indeed, it is easy to see, by induction on /c, that for any k <n and any particle 
j ^ + , the locations of j in id • Sj^ • • • Si^, and in Ti-Si,--- 5^^. are equal, 
since particle j does not distinguish between i and i + 1 when attempting a 
swap. For k = n, this shows that Y~^ is one of and TiX~^ . However, 
in Y, we must have Y~^{i -|- 1) < Y~^[i), since that is so in Tj and particles 
i and i + 1 cannot be swapped again when applying , . . . , Si„ . 

The lemma is now proved by induction. Assume it is true for some se- 
quence ii, . . . ,in and extend the sequence by adding io = i- We then have, 
applying the induction hypothesis for X (with X, Y as above) , that 

{id - SiQ - - - Si^) ^ = Y ^ = X - Sif,= id - Si^ ■ ■ ■ Sig, 

as claimed. □ 
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3. The infinite oriented swap process and exclusion processes. While 
our results are concerned with processes taking place within a finite interval 
[l,n], in order to prove them, it will be useful to see the interval as part 
of Z and the finite-interval process as a function of a process on the entire 
line. We now introduce this process, which will be called the infinite oriented 
swap process. 

In the infinite oriented swap process, particles with labels from Z occupy 
positions in Z, with each position containing exactly one particle at any 
given time, so the configuration space is Z^. The infinite oriented swap 
process is the Z^-valued continuous-time Markov process {Ct)t>o defined as 
follows. The initial state is the infinite identity configuration id [defined 
by id(fc) = k for all A; G Z]. For each edge {k, A; -|- 1), the particles in positions 
k and k + 1 "attempt to swap" at rate 1, succeeding if and only if they are 
in increasing order. 

The existence of this process is proved using a graphical representation. 
Specifically, for each k £Z, let Hi. be a Poisson process with density 1 on 
M"*" = [0,oo), where (n.k)k£Z form an independent family. For each k gZ, 
let Hk be the set of times at which a swap is attempted between positions 
k and A; -|- 1. Since almost surely, for all G Z and t > 0, we have that 
\Ilk n [0,t]| < oo, and, a.s. for all t > 0, we have lllj n [0,i]| = for some 
arbitrarily large (both positive and negative) values of j, it follows that the 
label of the particle in position k and the location of particle j at any time 
t are well-defined. 

As before, we call the configuration at time t, we call (^^{k) the location 
of particle k at time t and we call the function t C,^^{k) the trajectory of 
particle k. 

Note that for finite n, the oriented swap process {'rjt)t>o can be real- 
ized similarly to the infinite process by using only the Poisson processes 
\J^k}i<k<n-i s-iid ignoring all the others. More generally, for any (possibly 
infinite) interval / = [a, h] C Z, where a <h and a, & G ZU {itcxo}, we consider 
the oriented swap process on I defined by restricting configurations to / and 
applying only swaps coming from the Poisson processes {Hk}k&[a,h-i]- We 
denote this process by {Ct)t>o, with the conventions that if / = Z, then the 
superscript / is omitted, and that = c]^'""' = r/". 

The totally asymmetric simple exclusion process (TASEP) is a process on 
Z (or, for our purposes in some cases, a subinterval of Z) with just two kinds 
of particles (called a "particle" and a "hole"), where, with exponential rate 
1, each particle tries to jump one step to the right (leaving a hole behind), 
succeeding if the place to its right contains a hole. On a finite interval, 
particles cannot leave the interval from the right or enter it from the left, and 
the process will support only a finite number of moves before the particles 
are stuck at the right side of the interval. 
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It is well-known and easy to see that a TASEP on / = [a,b] can be con- 
structed starting from an arbitrary initial configuration using a family of 
independent Poisson processes {n.k} ke[a,b~i] of attempted jump times, sim- 
ilarly to the construction of the oriented swap processes above. We shall be 
interested in TASEPs with a particular class of step functions as the initial 
conditions. For any interval I = [a,b] Q7j and /c G Z, let (t'f '^)t>o denote the 
TASEP on the interval / with initial condition 

J^o^ix) = l{x<fc} 

constructed from the same infinite family of Poisson processes {Ilk}k^z that 
was used in the construction of the infinite oriented swap process above. As 
before, if / = Z, we may omit / from the superscript and denote the process 
simply as f^, and if / = we denote the corresponding process by t'f'". 

For any /c G Z and a configuration p G Z^, we define T^p G {0, 1}^ by 

{Tkp){x) = l{p(x)<k}, 
that is, the composition of the characteristic function of (— oo, k] with p. 

Lemma 3.1. With the above construction of the processes, a.s. for any 
subinterval / C Z and t> we have 



k,I 



Proof. The identity holds initially by definition and is preserved by 
any attempted swap. □ 

Our next goal is to describe the relations between (I'f '^)t>o and {i't)t>Q = 
{v^''^)t>o- Informally, if / C J C Z are two subintervals of Z (possibly infi- 
nite), then the randomness involved in (ff '^)t>o is a subset of the randomness 
involved in (f^''^)t>o. It turns out that I'f'^ is a function of f^'"^, not just as 
a process, but also at any fixed time t. 

For what follows, it will be more convenient to work exclusively with 
configurations in {0, 1}^ having a rightmost particle. Let C {0, 1}^ be 
the space of such configurations. If p G {0, 1}^, extend p to {0, 1}^ by setting 
p{x) = for X £ 7j\I. Note that Qq is a.s. invariant under the Markovian 
dynamics of the TASEP, so, because of our choice of initial conditions, we 
get that a.s. f^'^ G fio for all t > 0, A; G Z and / C Z. 

Define the operators Rk, and Jm (where m G Z) on as follows. We 
call Rk the cut-off operator. For a configuration p, Rkip) keeps only the k 
rightmost particles of p. Formally, 



{Rkp){x) 




if p{y) < k, 

y>x 

otherwise. 



12 



O. ANGEL, A. HOLROYD AND D. ROMIK 



The push-back operator Bn pushes all particles back into (— oo,n] and pre- 
serves the exclusion. moves the jth rightmost particle, if it is in location 
X, to location x A (n + 1 — j). Formally, 



{BnP){x) 



0, x> n, 

1, if x < n and ^ p{y) > n — x, 

y>x 

p{x), if X <n and ^ p{y) <n — x. 

y>x 



The jump operator Jm is analogous to the sorting operator Si from the 
previous section and tries to make a particle at m jump to m + 1 if there is 
a particle at m and no particle at m + 1. Formally, 



JmP- 



p- Tm, if p{m) = 1 and p{m + 1) = 0, 
p, otherwise. 



where p ■ denotes p with the values at m and m + 1 transposed. 

One more notion that will prove useful is the queue-length function. For 
any x € Z and a configuration p gQq, S{p, x) gives the number of particles 
to the right of x: 

Sip,x) = ^p{y). 

y>x 

Clearly, a configuration p gQq is completely determined by •). It is easy 
to see that in terms of the S'-function, the operators Rk and Bn take the 
following form: 

(4) S{RkP,x) = S{p,x)Ak, 

(5) S{BnP,x) = S{p,x)A{n-x) + 

(where, for it € M, we use the notation = uV 0). It follows that Bn and 
commute. We will also need the following result: 

Lemma 3.2. (i) // 1 < m < n - 1, then 

BnRk'Jm — J-mBnRk- 

(ii) If m>n, then 

BnRkJrn = BnRk- 

(iii) If m<0 and if p^VLq has its k rightmost particles in [1, oo), then 

BnRkJrnP = BnRkP- 
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Proof. We first prove (ii). i?„ and Rk commute, so it suffices to sliow 
that BnJm = Bn wlienever m>n. Tliis is true since Bnp depends only on 
S{p, k) for k <n, wliereas Jm does not change any of those numbers. 

Claim (iii) is similarly easy. RkP depends only on the positions of the k 
rightmost particles in p. Thus, if there are at least k particles in [l,oo) and 
m < 0, then does not affect any of them and RkJmP = RkP- Claim (iii) 
follows. 

To prove (i), we show that for m G [l,n — 1], Jm commutes with both Bn 
and Rk- First, we prove that RkJm = JmRk for any m G Z. Imagine that the 
k rightmost particles in a configuration p £ are colored red and that all 
other particles are colored blue. Then the statement that R^JmP = JmRkP 
simply says that making a particle at m, if there is one there, try to jump 
(note that its color will be preserved whether it jumps or not) and then 
deleting all blue particles is the same as first deleting all blue particles and 
then trying a jump at m. If p has no particle at m, then this is clearly true, 
and if it has a particle, then the statement can be false only if the particle is 
a red particle which has a blue particle to its right. But this cannot happen 
since all blue particles are to the left of all red particles. 

Next, we prove that if m < n, then BnJm = JmBn- If, in a configuration 
p, there is no particle at m, then leaves p unchanged. Then, if, in BnP, 
there is also no particle at m, we are done; otherwise, there is a particle that 
was pushed there by other particles from the right, so, in particular, BnP 
also has a particle at m + 1 (here, we use the fact that m <n), and Jm also 
leaves BnP unchanged, so we are done. 

Alternatively, assume that p has a particle at m. Let j denote the ranking 
of that particle (in terms of right-to-left order of appearance) and let m' 
denote the location of the jth. rightmost particle in BnP (we think of it as the 
"same" particle after the push-back operation). If m' < m, then the particle 
was pushed, so m' = n + 1 — j. Applying Jm will leave BnP unchanged since, 
in BnP, there are particles at all locations between m' and n, including m + 1. 
In the other direction, applying Jm to p might make the jth particle jump 
to m + 1, but, after applying Bn, it will again be pushed to m' = n + l — j. 
So, in this case, too, we have shown that BnJmP = BnP = JmBnP (note 
that all other particles are unaffected by Jm)- It remains to deal with the 
case m' = m: in this case, if m = n + 1 — j, then whether applying Jm to 
p produces a jump or not, the jth particle will be pushed back by Bn to 
n + 1 — j, with a particle in front of it blocking a jump (since m <n). And, 
if m < n + 1 — j, then the jth particle does not get pushed back, so, if p 
has a particle at m + 1, that particle will still be there after applying Bn 
and, in both cases, a jump will not occur; otherwise, if there is no particle 
at m + 1, after applying Bn, the (j — l)th rightmost particle will be pushed 
to n + 2 — j > m + 1, so, after applying Bn, there will still not be a particle 
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at m + 1 and a jump will occur for both orders of applying the operations. 
This completes the proof that BnJm. = JmBn and hence the proof of (i). □ 

Lemma 3.3. Let k,n G Z. Almost surely, for any t>0, we have 
(6) i/f'" = S„i2fci/f. 

Proof. The equality (6) is satisfied at time t = 0. Attempted jumps 
outside of [l,n — 1] have no effect on u^'"" and, by parts (ii) and (iii) of 
Lemma 3.2, these have no effect on BnRk^t either. Jumps inside [l,re — 1] 
a.s. occur at a discrete set of times and, by Lemma 3.2(i), whenever a jump 
is attempted between positions (m,m + 1), if (6) was satisfied before the 
attempted jump, it will remain true after it. □ 

4. Hydrodynamic limits. The TASEP {i^t)t>Q with initial condition Uq = 

l{-oo,o] ti^s been studied in great depth and is the simplest case of a shock 
in a TASEP. Lemma 3.3 allows us to tap into this knowledge. We use the 
following fundamental result of Rost ([13], Theorem 1; see also [10], Chapter 
VIII, Section 5): 

Theorem 4.1 (Rost [13]). For any — oo <u <v < oo, we have that a.s. 
lim - u^{j)= h{x)dx, 



ut<j<vt 



where 



{1, x<-l, 
1 — X 
-l<x<l, 
0, x>l. 

The following immediate corollary is an equivalent formulation of Post's 
result which can be interpreted as saying that the family of functions (x — > 
h{x/s))s>o is the limiting time-evolution, or hydrodynamic limit, of the den- 
sity profile of the process {v^)t>o, when one scales both the time- and space- 
axes by a parameter n that goes to infinity. 

Corollary 4.2. For any —oo < u < v < oo and s > 0, we have a.s. 

nu<j<nv 

With this preparation, we can now formulate and prove an analogous 
hydrodynamic limit theorem for the TASEP z^f where y ■ n for some 
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fixed < y < 1. In terms of tlie original oriented swap process {rj^)t^ this 
describes tlie limiting flow of particles with label < /c, which start out in 
positions [1, A;], to their final positions at [n + 1 — as a function of the 
scaled time parameter. This result is essentially an encoded form of Theorem 
1.2 and we will derive Theorem 1.2 from it later. 



Theorem 4.3. Fix 0<y<l, s>0 and Q <u <v <l. Let k = k{n) he 
a sequence of integers such that k/n ^ y as n ^ oo. Then, the number of 
particles in the oriented swap process {rfl)t with index < k that are in the 
interval {nu,nv) at time us satisfies 

(7) i#{l < i < A: : (O-^j) G {nu, nv)] £ fis, x, y) dx, 

where f is the function of s,x and y defined by 

^ + y-^ {y-s)y L~{s)<x<{y + s)AL+{s), 



(8) fis,x,y) 



2s ' 

1, {) < X <y — s, 



0, y + s < a; < 1, 

1, s> I- y and {I - y)\/ L^{s) < X <l, 
0, s> y and d < X < {1 — y) /\Ly {s), 

where L^{s) are as defined in Theorem 1.1 in the Introduction. 

Figure 6 shows two convenient ways of visualizing the limiting density pro- 
file /(s, X, y) for fixed y (in this case, y = 0.3): in Figure 6(a), we see /(s, x, y) 
in the (s,x)-plane; Figure 6(b) shows a succession of plots of f{s,x,y) as a 
function of x, for several increasing values of s. If we think of particles with 
label < A: in the oriented swap process as "red" particles and particles with 
label > A; as "green," this illustrates how the red particles advance into the 
green zone, eventually displacing all green particles in positions [n + 1 — A;, n] . 

Proof of Theorem 4.3. It is clearly enough to prove (7) for v = oo. 
li k/n^ y as oo, then, by the definitions, 

-#{1 < 3 < k:{nlsr\j) e {nu,^)} = -5(i.^r,nn). 
n n 

By Lemma 3.3, this is a.s. equal to 

i[5'(t'^„,nn) A A; A (n - |nnj)+]. 
n 

By Corollary 4.2, together with translation-invariance of the TASEP dy- 
namics, as n ^ oo, this last quantity converges in probability to 



(9) F{s,u,y):=J h(^^^-^^ dx A y A {I - u) 



+ 
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(a) (b) 

Fig. 6. (a) Schematic representation of f(s,x,y) in the {s,x)-plane for j/ = 0.3; (b) 
several time slices of f{s,x,y) (shown as a function of x) for y = 0.3. 



It now remains to show that 

(10) F{s,u,y)= f{s,x,y)dx, 

J u 

where f{s,x,y) is defined in (8). Define the functions 

A~{s) = inf |n : £ ^(^) > y}, 

A+(s)=inf|M: h(^^^^dx>{l-u) + y 

We can evaluate A~ and A+ exphcitly, as shown in Figure 7, to give 

JO, 0<s<y, 



1, 0<s<l-y, 
LUs), s>l-y 



[e.g., from Figure 7, it is easy to see that A~(s) for s > y is the solution of 
the equation (y + s — A)^/2s = y; we omit the detailed verification of the 
above formulae]. 

Now, check that A^(s) < A+(s) if < s < and A~(s) > A+(s) if s > 7j^. 
Equipped with this information, it is easy to write the following explicit 
formulas for F{s,u,y). Assuming that y < 1/2, we get 

y — u, < n < y — s, 

F(^,n,y) = <i i^ + y-^)\ y-s<u<y + s, (ifO<s<y), 
4s 

0, y + s<u<l 
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Fig. 7. The functions Ay (s), A+(s). 



0<u<Ly{s), 

Fis,u,y) = \ ll±|^, L~{s)<u<y + s, (ify<s<l-y), 
0, y + s<u<l 

y, o<u<Ly{s), 

F{s,u,y) = \ (^ + ^~")' , Lyis)<u<L~{s), {if I - y < s <jy), 
l-u, Ly{s)<u<l 

J?/ \ 0<U<l-y, /-r ^ ^ 

^(^'^'^) = ll-n, l-y<u<l (^f^>^^)- 

On the other hand, referring to Figure 6 for convenience, one may rewrite 
f{s,x,y) defined in (8) more explicitly (again in the case y < 1/2) as 

1, < X < y — s, 

f{s,x,y) = { — , y-s<x<y + s, (if < s < y), 

0, y + s<a;<l 

0, 0<x<Ly{s), 

f{s,x,y) = \ ^ + y-^ ^ L-is)<x<y + s, (ify<s<l-y), 

0, y + s<a;<l 

0, 0<x<Ly{s), 
f{s,x,y) = \ ^^y^ ^ , Ly{s)<x<Lyis), (if 1 - y < s < 7^), 

1, Ly{s)<X<l 

^(^'"'^) = ll, l-y<x<l (^f^>^^)- 

Comparing the two sets of formulae, it is clear that (10) holds. The case 
l/2<y<l is dealt with similarly and is omitted. □ 
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Proof of Theorem 1.2. To prove (2), it is enough to prove that for 
each s > and < x, ?/ < 1, we have that 

M^([0,x] X [0,y])-^Ks{[0,x] X [0,y]). 

However, looking at the definitions and using Theorem 4.3, we see that 

/i:([0,x] X [0,y]) = -#{l<j<ny:{r,:i,r\j)e[0,nx]} 

n 

r 

f{s,u,y) du. 

It 'l_XJ 

It remains to verify that 

f{s,u,y)du = Ks{[0,x] X [0,y]) 

'0 

or, equivalently, that 

d 

(11) f{s,x,y) = —Ks{[0,x] X [0,y]) =:g{s,x,y). 

Note that, for s > 2, Kg is simply the arc-length measure on the segment 
{{x,y) G [0, 1]^ :y = 1 — x}, normalized to be a probability measure, so the 
identity is trivial to check. There are two other cases, < s < 1 and 1 < s < 2, 
and in each of these cases, it is easy to compute g{s,x,y) by dividing the 
unit square into the various possibilities for {x,y). Figures 8(a) and (b) 
show the result (see also Section 6, where we give an alternative geometric 
description of Ks). Comparing this to the explicit formulas for f{s, x, y) in the 
case y < 1/2, one verifies that (11) holds in this case. The case 1/2 < y < 1 
may be dealt with similarly or can be deduced from the symmetry of the 
process with respect to reversing left and right, and replacing particle k by 
n + l-k. □ 



5. Finishing times. 

Proof of Theorem 1.5. To analyze the time /3"(A:) at which particle 
k completes its movement, note that, in the oriented swap process {r]^)t, 
particle k stops moving as soon as the particle is in position n + 1 — k 
and particles {l,...,k — 1} are to its right. Define the events (depending 
implicitly on n) 

Al = {S{u^''',n-k)=k} 

= {particles {1, . . . , A;} are in positions {n — A; + 1, . . . , n}}. 
The key to the result is the identity 

(12) {p''{k)<t} = A'^nA'^-\ 
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(a) (b) 

Fig. 8. Values of g{s,x,y) in the {x,y) plane for (a) < s < 1 and (b) 1 < s < 2 (there is 
no need to treat the boundary case s = 1 because of continuity in s ). To compute g{s, x, y) , 
integrate Kg over the shaded thin strip of infinitesimal width. 



By Lemma 3.3, we can write 

A'l = {S{BnRkiytn-k) = k} 

(13) = {S{v^,n -k) AkA{n-{n- k))+ = k} 

= {S{u^,n-k) > k}. 

Essentially, the fact that S{u^,-) has a deterministic scaling limit implies 
that has a sharp threshold (that is, its probability increases from e to 
1 — £ over an asymptotically small time interval). Thus, the finishing times 
are concentrated. 

By symmetry, it suffices to prove that, with high probability, P^(k) is 
close to "fk/n foi^ all 1 < A; < n/2. With this in mind, we start by showing this 
for a single particle. Let k = k{n) be such that k/n ^ y G [0, 1/2] as n ^ oo. 
From Theorem 4.3, we have, for any fixed s, that 

— S(u^^,n—k) — - — > [ hi -] dx 

n ^ n-oo \ s J 

(14) 

r 0, < s < 1 - 2y, 

I 4^ ' 

Thus, and A^~^ hold with probability tending to 1 as soon as > 
y, which simplifies to s > 1 + 2^/y{l — y) = 7^. Conversely, if s < 7^ and 

y / 0, then P(^^J and P(A^;^) tend to 0. This implies that /3"(A;)/n ^ 7^ 
for any y G (0, l). 

To get simultaneous convergence for all particles, we will use the easy 
facts that S{v'^g,x) is increasing in k and s, and decreasing in x. Since, for 
any fixed s,y, we have convergence in probability of n~^S{i'^g,n — k), this 
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monotonicity implies that n~^S{u^g,n — k) converges in probability to the 
right-hand side of (14), uniformly in k and s < 2. [Proximity to the limit at 
a finite set of implies proximity at all intermediate points.] 

What this implies for the finishing times is that for any e > 0, the scaled 
finishing times /3"'(fc)/n for particles /c G { \_en\ , . . . , [n/2j } are (with asymp- 
totically high probability) uniformly close to the given limit 7fc/.„, that is, 

P 

max |/3"'(A;)/n — 7i,/„ I >0. 

For A; G {1, . . . , [enj}, however, we only get an upper bound on the finishing 
times, namely, that 

max (/J"(A;)/n-7,/„)+^0. 

Note that 7^/^ is close to 70 = 1 for such small fc, so we know that the 
particles with small labels must finish shortly after time s = 1 , but not that 
they cannot finish much sooner. The reason we do not yet get a lower bound 
when k <^n\s that for such /c, the scaling limits of k and S^u^f^-^^^-^ ,n — k) are 

k 



both 0. Thus, the scaling limits are not enough to deduce that S{i'^,^ 



(l-e) 



n - 



k) < k and hence that the event A^g occurred. To complete the proof, it 
suffices to note that particle 1 performs a random walk, moving only to 
the right at random times with rate 1, and, therefore, at time (1 — 2e)n, 
it is, with high probability, still at a distance of at least en from position 
n. However, this implies that particles {1, . . . , [enJ} are also not finished by 
time (1 — e)n. □ 

Proof of Theorem 1.6. For a configuration p£ Qq, let vrp(j) denote 
the position of the jth rightmost particle in p. By (12) and (13), we see that 

/3"(A;) =y"(A:) vy"(A;-l), 

where 

y"(j)=inf{t>0:5(z./,n-j)>j} 

= inf{t>0:7r^.(j)=n + l-j}. 

If we were interested in V"'{k), it would follow immediately from a theorem 
of Johansson [8], Theorem 1.6, that V"'{k) converges in distribution after 
scaling to the Tracy- Widom distribution. Because our random time is a 
maximum of V"'{k) and V"'{k — 1), we need to show that these two times 
cannot be very far apart. Define random times Wi = Wi and W2 = W2 by 

H^i = inf{t > : vr u-Ak - I) = n + 1 - k}, 
W2 = inf{t > : vr fc-i (A;) = n + 1 - A;} 
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and observe that we have the bounds 

Wi < (3'^{k) = Vik) V T/"(A; - 1) < W2. 
The lower bound follows from the fact that 

7T^k{k)<Tr^k-i{k-l), 

which implies that V"'{k) > Wi. Similarly, the upper bound follows from the 
inequalities 

< ~ 1) 7r^fc-i(A:) < 7r^fe(/c + 1) < 7r^fe(A:), 

which imply, respectively, that V^{k — 1) < W2 and V''^{k) < W2- Now, use 
Johansson's theorem ([8], Theorem 1.6, see also [8], Corollary 1.7) for Wi 
and W2, to get that 

(15) —^57:5 ^> Ftw as n — > 00 (i = 1, 2) 

(we use the well-known equivalence between the TASEP with Rost's step 
initial conditions l(_oo,o] s-'^d directed last-passage percolation in N^; in the 
notation of the paper [8], we have N = yn and 7 = (1 — y)/y). Since f3"{k) is 
bounded between two random variables having the same distributional limit 
(with the same scaling), it must also converge in distribution to .Ftw with 
the same scaling. □ 

6. The inversion number. In this section, we prove Theorem 1.4. By 
Theorem 1.5, it is enough to consider s < 2 since, for s > 2, we have, with 
high probability, that 77"^ is the reverse permutation and, in particular, 

mv(CJ = (2)- 

Let r] £ Sn be any permutation with normalized empirical measure fi de- 
fined as in (1). Thus, /i is supported within the square [0,1]^. For points 
z = (x, y) and z' = (x' , y') in [0, 1]^, we write {z \ z'} if {x < x' and y > y'}. 
The basic observation is that we can express the inversion number inv(r7) in 
terms of the measure fi, as follows. If we sample independent random points 
z,z' £ [0, 1]^ with distribution fi, then 

inv [t]) = IJ,{z\z') = [ [ l[z\ z']fi{dz')fi{dz). 



'[0,1]2 J[0,l]2 

The integral above, as a function of the measure /i, is continuous at all 
measures that assign measure to all horizontal and vertical lines (since, 
for such fi, the product measure will assign measure to the boundary 
of the set {z \ z'} C [0, 1]^). This includes, in particular, the measures Kg. 
By Theorem 1.2, it follows that 

invivls)-^Iis):= f / 2.1[z\z']Ksidz')^sidz). 

\2y n^oo J[0,l]2 7[0,1]2 
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Thus, calculating the limit of the scaled number of inversions function is 
reduced to evaluating the Kg (8) Ks measure of a certain set in [0,1]^. This 
integration is tricky, so we sketch an argument below. 

We first give a geometric description of Ks which will shed some light 
on subsequent formulae. To sample a point {x,y) from Kg, first choose y 
uniformly in [0, 1] and x uniformly in [y — s,y + s] [so that (x, y) is uniform 
in a parallelogram of area 2s]. Next, evaluate L~{s) and L^{s). If x is outside 
the interval replace it by whichever of L~ is closer to x. When 

s > 1, if L~ > L'^ , then replace x hy 1 — y (see Figure 9). The equivalence of 
this description to the original definition of Kg can be verified with a simple 
computation that we omit (it also has a more conceptual explanation in 
terms of Theorems 1.2 and 4.3 and the combinatorial operations Bn and 
Rk] e.g., look again at Figure 7). 

Claim 6.1. For a point z = (x, y) in the support of Ks, if z' has law Kg, 
then 



Proof. The second-case is trivial. The first case is equivalent to showing 
that the part of the parallelogram with z' satisfying z\ z' has area {y — 
X + s)^/2. This region is a fairly simple polygon [see Figure 10(a)] and its 
area is easily computed. 

When s < 1, there is a nice geometric argument. Since Kg is symmetric 
with respect to the refiection {x,y) — > {y,x) (check directly, or deduce it 
from Theorems 1.2 and 1.3), the symmetric description of how to sample 
z' by choosing x' first and then y' , and shifting y' vertically if the resulting 
point is outside the support of Kg is equally valid. Therefore, we can reflect 
the curved triangle ABF of Figure 9 and deduce that Kg{z' : z \ z'), which 
is equal to the area of the shaded region in Figure 10(a), is also the area of 
the shaded region in Figure 10(b), which is simply a right triangle. 




s<l orye[Wg-,Wg+], 
s>l andyi [W' ,W+]. 



D 



C 



D 



C 



A- 




A- 



B 



B 



Fig. 9. The measure KsI a point is selected uniformly in the parallelogram ABCD . If it 
is outside the curved polygon BEDF where Ks is supported, then it is shifted horizontally 
into the support. A similar description also applies for s > 1. 
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(a) (b) 

Fig. 10. Proof of the claim when s < 1. We will have z' \, z if the point from which z 
is obtained by shifting the x-coordinate lies in the shaded region in (a). Symmetry of Ha 
implies that the areas of the shaded regions in (a) and (b) are equal. 



For the case s > 1, we may either deduce the claim from the case s < 
1 together with analyticity of the areas involved, or prove it directly by 
computing the area of the hexagon, as in Figure 10(a). □ 

We now compute I{s), starting with the case s < 1. From the claim, it 
follows that, in this case, 

I[s) = / 2 Ks[dz). 

J[o,i]^ 4s 

Here, the contribution from z on the lower arc [with x = L~{s)] is 



2 

dy = — . 

Jy-s 4s 2s 5 

The upper arc gives an equal contribution. 

To integrate (y — x + s)^/4s in the interior, note that, for any s < 1/2, 
the integral of a function G in the internal region is given by 

(16) / / Gdxdy+ / / Gdxdy+ / / Gdxdy. 

Jo J Js Jy—s Jl~sJy~s 

For 1/2 < s < 1, it seems, at first glance, that a different formula is required, 
but, in fact, it is easy to see that (16) is still correct due to cancellation. 
Applying this with G = {y — x + s)"^ / (8s^) and adding the arc contributions 
gives the answer /(s) = |s — ^s^ for s < 1. 

The case s > 1 is trickier. The contributions from z with y < W~ or 
y > are simple enough. The integrals for the arcs and interior are also 
not difficult to evaluate, but, due to having as endpoints, give less 
elegant expressions. Overall, the different parts shown in Figure 11 evaluate 
to 

A + A' = 1- X, 

B = B' = X + -l-s-^/'^{{l - Xf^ - (1 + Xf/^), 
5\/2 
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Fig. 11. Contribution to I{s) for s > 1. A, A': diagonal segments; B,B' : arcs; C: 
lens-shaped region between the arcs. 

c = (-1 - - - xfi^ - (1 + xfi^) 

where we denote 

X = ^s{2-s), Y = vT^^ - VlTX, 

z = VT^ + ViTx. 

It seems impossible to simplify this any further until one realizes that Y = 
— ^2(2 — s) and Z = \f2s. Then, straightforward manipulation of I(s) = 
A + yl' + B + B' + C gives the theorem. 

7. Particle trajectories and second-class particles. Our goal in this sec- 
tion is to prove Theorem 1.1. To this end, note that one can identify the 
location of particle k in the oriented swap process {r\i \ at a given time i as 
the unique place where the two configurations and i/f differ. Because 
of Lemma 3.3, this is the same as the place where the two configurations 
BnRki^f and BnRk-i^t'^ differ. 

Note that f^, z^f ~^ also differ in a unique place. This is true at time t = 0, 
and it is easy to see that it remains true whenever one of the jump operators 
Jm is applied. In the context of the TASEP, the place where v^, i^t~^ differ is 
referred to as a second-class particle since it behaves in effect like a particle 
which has priority over holes (i.e., can swap with a hole to its right and 
does so with exponential rate 1), but over which the "first-class" particles 
initially at (— oo. A; — 1] have priority (to see this, try to imagine the effect 
that the operator has on the pair of configurations). 

As before, let denote the space of configurations in {0, 1}^ having a 
rightmost particle and, if a configuration is in {0, 1}^ for some subinterval 
/ C Z, consider it to be in {0, 1}^ by extending its definition to be outside 
I. We call two configurations p, p' G ^Iq compatible if p and p' differ in exactly 
one position where = p'{x) < p{x) = 1, and we denote the position where 
they differ by ^p,p'- For p £ ^Iq, denote the position of the jth rightmost 
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particle of p by vrp(j) and denote by 9p{n,j) the position of the jth rightmost 
hole of pV l[„,+i,oo) (i'^ words, the jth rightmost hole among the holes that 
are in ( — oo , n] ; define it as n + 1 if j = ) . The following lemma elucidates the 
effect that the transformations Bn and have on the second-class particle. 

Lemma 7.1. Let p,p' G f^o be compatible configurations. Then the pair 
RkP, Rk-ip' is compatible, with its second-class particle in position 

(17) ^RkpA_,p' = ^p,p' V TTpik). 

If p has infinitely many holes in (— oo,0], then the pair Bnp,Bnp' is also 
compatible, with a second-class particle in 

(18) ^B^p,B^p' = ^p,p' A ep{n,S{p, n)). 

Consequently, for k <n and any compatible p,p' E Qq, the pair of configu- 
rations BnRkP, BnRk-ip' is compatible with a second-class particle in 

(19) ^B„Rkp,B,,Rk.^p' = T.pyV TTp{k) AeR^p{n,S{Rkp,n)). 

Proof. Note that (19) follows from (17) and (18) [since RkP satisfies 
the assumption of (18)]. For (17), note that if ^p,p' > T^pik), then, when 
applying Rf^ to p and also when applying Rk-i to p' , all particles to the left 
of TTp{k) = 'Kpi{k — 1) will be deleted, so the transformed configurations still 
differ at Sp p/. If Sp p/ < T^p^k), then, in RkP, all particles to the left of TTp^k) 
are deleted and in R^-ip' , all particles to the left of iTp'{k — 1) = 7rp(fe — 1) 
are deleted, so the two transformed configurations now differ at T^p{k). This 
proves (17). 

To prove (18), first note that an equivalent description of the transforma- 
tion Bn acting on p is that it takes the S{p, n) particles to the right of n and 
places them in the 5(p, n) rightmost holes to the left of n + 1. This is true for 
configurations with sufficiently many holes to the left of n + 1, which is why, 
for convenience, we assumed that there are infinitely many such holes in p. 
Similarly, Bn acting on p' has the effect of filling in the S{p',n) rightmost 
holes in p' to the left of n + 1. Now, assume that Sp^p/ < 9p{n, S{p,n)). In 
particular, it follows that ^p,p' < n + 1, so S{p,n) = S{p',n). Because of the 
assumption, the S{p,n) rightmost holes to the left of n + 1 are the same for 
p and p' , and they both get filled under the operation Bn- So, BnP and Bnp' 
still differ only at position 5^p,p' and (18) is correct. The other possibility is 
that Sp^p/ > 9p{n, S{p,n)) [we cannot have equality since p has a particle at 
position Spy and a hole at position 6p{n, S{p,n))]. For this case, consider 
two further sub-cases: first, if ^p,p' < n, then S{p, n) = S{p' , n), and then the 
holes in p in positions 6p{n,l),9p{n,2), . . . ,6p{n, S{p,n)) get filled in BnP, 
whereas for p' , the holes in positions 6p{n, 1), . . . , Op{n,S{n, p) — 1), Ep^p/ get 
filled when applying Bn- It follows that BnP and Bnp' differ in position 
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6p{n, S{p,n)), as claimed. Finally, if Sp^p/ > n, then S{p',n) = S{p,n) — 1, 
and then, in p, the holes in positions 9p{n,l), . . . ,0p{n, S{p,n)) get filled, 
but in p' , the holes in positions Op{n, 1), . . . , 6p{n, S{p, n) — 1) get filled, so, 
again, BnP and Bnp' differ in position 9p{n, S{n, p)). □ 

Denote by Xt the location of the second-class particle in the TASEP pair 
i^^,!/^^. Guiol and Mountford [11] proved that Xt/t converges a.s. as t — > oo 
to a random variable distributed uniformly in [—1,1]. The following result 
is an easy consequence of their result. A similar, but slightly weaker, result 
was proven earlier by Ferrari and Kipnis [7]: there, the limit is with respect 
to the (weaker) topology of pointwise convergence of functions. 

Theorem 7.2 (Guiol and Mountford [11]). Let U denote a random vari- 
able distributed uniformly on [—1,1]. Let X^{s) = n~^Xns denote the tra- 
jectory of the second-class particle when space and time are scaled by n, 
considered as a random function in the function space IR['''°°) . Then, 

X^(s)^U-s. 

Here, denotes convergence in distribution with respect to the uniform 

topology on functions [0, cxo) M. 

Proof of Theorem 1.1. Let k = k{n) be such that k/n^ y € [0, 1] 
when n — > oo. For s > ^y, we already know, from Theorem 1.5, that with 
asymptotically high probability, TJ}{s) = — > 1 — y = (/)j,(s), so it will be 
enough to prove the claimed convergence on the space of functions M^''^^']. 
From Lemmas 3.3 and 7.1, we have that 

Tu (s)= X] k.n fc — l.n= d k D D A; — 1 

= V7r^fe^(A;) A6'^^^fc^(n,S'(i/^,,n) A/c)]. 

Therefore, we need to understand each of the components on the right- 
hand side. Note that, by translation, S ^ k-i is equal in distribution to 

k -\- S o ,,-1 = k-\- X'^(s). Therefore, by Theorem 7.2, we know that 
We now claim that 

(20) n-\,.Jk) ^ A- is), 

n~^S{v^g,n) — — > ■ipy{s) := J dx 
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(21) 



(see Fig. 7) 




0<s< 
s>l-y, 



(22) ej,^^,Jn,Siul,n)Ak) ^ 
If we prove these claims, it will follow that 

(23) ^ (y + Us) V Ay{s) A A+(s) = <A,(s) {s < -fy), 

which is the claim of the theorem. Relation (21) follows immediately from 
Corollary 4.2. The other two relations, (20) and (22), are also relatively 
straightforward consequences of Corollary 4.2. To prove (20), note that 
Tr,,k (k) satisfies 

ns 

S{'^ns,^u^Jk)) = k-l 

[indeed, it is the minimal x with S{v^^g,x) = k — 1]. On the other hand, if 
y > 0, then, by Corollary 4.2, for any e > 0, we have that, as n — > oo, 

n-i5(i/i,(A-(s) + e)n) 

(24) 

h( ] dx <y — 6 and 

JA-{s)+e \ s J 



. nsi \ y 

25 

roo 



h(- — -] dx>y + 5 

Ayisye \ S J 

for some 5 = S{y,£) > that depends on e. By monotonicity of S{p,x) in x, 
it follows that 



F[A-{s) - e < n-'TT^M < A^is) + e] ^ 1 



as n ^ oo, 



which implies (20) since e was arbitrary. 

In the extremal case y = (which, by symmetry, also implies the case 
y = 1), (25) still holds, giving the lower bound for n^k (k), and (24) does 
not, but we instead use the easy fact that 



(s) as n ^ oo 



(since the rightmost particle performs a random walk jumping to the right 
at rate 1) to get the upper bound, so (20) still holds. 

Next, to prove (22), we first simplify the left-hand side by noting the 
easily checked fact that ipy{s) < y for 0<s<l-|-y-|- 2y^ and, in particular, 
for < s < 7y since 7y < 1 + y + 2y^. Therefore, we have 

hs{ut,n)Ak)-^^y{s). 



fl n— >oo 
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Define 

fO, x<A-{s), 

By Corollary 4.2 and (20), Zy^s{x) is the limiting density profile of the pro- 
cess RkU^g, in the sense that, for all n € M, we have 

n^^S{RkVns^jn) ^ > / Zy^s{x)dx. 
n~*oo 7„ 

It follows, by an argument similar to the one used to prove (20) above, that 
n~^^R^^k^{n, S{i^^„n) A k) 

^^^> infjn: J {1 - Zy^s{x)) dx > ipyis)^ . 

For s < 7jy, this last function is easily seen (refer again to Figure 7) to equal 
Ay{s). This completes the proof of (22) and therefore also of (23). □ 

8. Additional comments and open problems. 

1. Uniformly random sorting networks. The uniform sorting network model 
of [2] exhibits behavior similar to the oriented swap process. A key problem 
is to prove the conjectures in [2] that are the analogs of Theorems 1.1 and 
1.2. One possible approach would be to try to relate uniform sorting net- 
works to some variant of the random swap process which can be analyzed 
using the exclusion process techniques developed in this paper. 

2. Limiting distribution of the absorbing time. Theorem 1.6 gives the lim- 
iting distribution of the fluctuations of the finishing times of individual par- 
ticles. However, the relation between finishing times of different particles is 
more delicate and requires knowledge about the joint distribution of last- 
passage percolation times. An interesting open problem would be to find 
sequences of scaling constants {an)^=i, {bn)^=i and a distribution function 
F such that the absorbing time /?" of the oriented swap process satisfies the 
convergence in distribution 

On(/3" — 2n) — bn =^ F as n — > oo. 

3. Partially asymmetric swap processes. The asymmetric exclusion pro- 
cess (ASEP) is defined similarly to the TASEP, with the difference that 
particles can jump in either direction, where jumps to the right happen at 
rate a and to the left at rate a' . When q' = 0, this is (a time-change of) the 
TASEP. When a > a' > 0, this is the partially asymmetric exclusion process 
(PASEP). 
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Some results for the TASEP have known analogs for the PASEP, while 
others do not. In particular, the following conjecture is needed to prove an 
analogue of Theorem 1.1 for the partially asymmetric swap process men- 
tioned in the Introduction: 

Conjecture 8.1. Theorem 7.2 holds for the partially asymmetric model 
with parameters a > a', with speed U uniform on [a' — a,a — a']. 
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